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ABSTRACT 


This report consists of two major parts. The first deals vith the 
development of formulas for computing the probability that a point taken 
from @ normal bivariate elliptical distribution with specified mean 
(X,Y ) end standard deviations (d, , Sy) shall fall within « circle of 


given radius (R) whose center is displaced a given distance (D) from the | 
center of the distribution. The second purt consists entirely of probebdil- 
ity tables. The entries (input parameters) to these tables are given in 
units standardized vith respect to 6, , mmely: RM, , B/6,, a . the 
spectra of values for these input parameters are as follovs: 

Oy MS, :0.0(0.2) 4.03 R/d,:0.0(0.1)%5 3 D/Sx: 0.00(00.05) 5.00. 

These tables will prove especially useful in dealing with problems involving 
accuracy studies of weapons systems and with other problems notably in 
meteorological studies. The events in many practical probability problems 
are best described by a normal bivariate elliptical distribution with 
unequal standard deviations. For example, one may be confronted vith the 
problem of evaluating the probability that a missile will hit a circle of a 
specified radius whose center (aim point) is displaced a given distance 
from the mean (of impact points) of a normal bivariate elliptical distribu- 
tion. In this example the impact points are governed by a noraal(Gaussian) 
bivariate elliptical density function; the mean of this distribution is not 
zero (1.@., the center of the distribution is not about the aim point). 


ii 
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2. SOOARY 
3. DIscUSsIG® 
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bd. Offset Circle Probabilities for the Normal 


Bivariate Elliptical Distribution 
c. Examples of the Use of the Tables 


TABLES 


Offset Circle Probabilities for Normal Biveriace 
Elliptical Distributions 


- 


6B wn 


GLOSSARY 


N[(=y); (2,9) 5 (6, ) Oy) ] : Two-dimensional density function of a 
. ‘normal (Gaussian) bivariate elliptical 
, distribution with zero correlation and 


| i 
(X%¥) : Coordinates of the mean of the distributica. 
(6, ,5y) 3 Standard deviaticas of the x and y variates. 
M(z%,y) : Mean ("aim dias”) of the distribution. 


5 : Any region in the xy plane. 


fe : Cireular region of radius R centered at origin (0,0) ("aim point”) 
and offset a distance D from the mean M (X,Y) of @ normal bivariate 
elliptical distribution. 


R : Radius of circle centered at aim point (0,0). 


D : Offset distance from center of circle to mean of the distribution 
(= R454) | 
F[(x,y)€S] : The probability thet « point (x,y) from a normal 
bivariate elliptical distribution N[(xy (2,95 


O,# Sy] with zero correlation falls vithin a 
region S in the xy plane. . 


P [cz,y) € £} : The probability that point (x, y) from a normal 
bivariate elliptical distribution N[(zx,y); (Z,7); 
= dy ] falls within a circular region of 


radius R lying in the plane along an axis of symmetry 
centered about the origin (0,0) and offset a distance 
D from the mean (X,Y) of the distribution.. 


r : Ratio of d, to 6, (=, /4,54, #0) 


RM, , Dh ,5,/5, + Dimensionless stendardized units vith respect to 


Input parameters (entries) for the probability 
tables appearing in this report. 


A: R/6, 
B: Diy 


GLOSSARY (Contd) 
a: ar pyetsurt =r’ 
do: rt 
o: (Rng) =R/26, =A /ar? 
, p¥262= 4)" 
: (72-1) cos*@ =bcos*@ 
:2r7p cos 8@=% cos 6 
:2r%6? = 20% 
: 2r-p=Vac 


(p,@) : Polar coordinates to which rectangular coordinates (x, y) are 
transformed. 


2a vy > @ 


exp( ) :e6 ) (e = Naperian base of natural logarithms) 
>_: Symbol for summation (F XE Rte In) 

x,y : . a =a é)= a(x, ) 
J zs) Jacobian of (x, y) with respect to (p,@)= Stee 
n! Factorial n (#1.2.3...n) 


C. _"t Nat)eee t+ n' 
~ i! U(n-t)! 


Ci): Binomial coefficients = nCirCyis 


0.0(0.2)4.0 : 0.0, 0.2, 0.4, 0.6, e +9 4.0 e 


TABLES OF OFFSET CIRCLE PROBABILITIZS FoR 
A NORMAL BIVARIATE ELLIPTICAL DISTRIBUTION 


1. INTRODUCTION 


Need for Offset Circle Probability Tables 


Apart from considerations of a purely theoretical nature, many 
practical (probability) problems present themselves with reference to events 
vest described by a bivariate normal distribution with unequal standard 
deviations (é.e. oy = Oy ) -» The standard deviations in the x and y 


coordinates are represented by dO, and Oy respectively. Such a distribution 
is known as a normal bivariate elliptical distribution. 


One may, for example, be confronted with the problem of evaluating the 
probability thet a missile will hit a circle of radius of a specified magni - 
tude whose center (aim point) is displaced a given distance from the mean 
(of impact points) of @ normal bivariate distribution with Oy # dy . 

Note that in this particular eyvample the impact points are governed by a 
normal (Gaussian) bivariate elliptical density function; the mean of this 
distribution is not zero (i.e. the center of the distribution is not about 


the aim point). 


These tables will be especially useful in providing ancwers to such 
problems as weli’as to kindred problems dealing with accuracy studies of 
weapons systems, their capability of delivering missiles with a specified 
accuracy, etc. (e.g. accuracy studies of the aim point and impact point). 


There sre, indeed, other problems, nctably in meteorological investiga- 
tions, where the solutions to problems are best approximated and dezvribed by 
normal bivariate elliptical distributions. 


2. SUMMARY 

This report consists of two major parts. The first part deals with the 
development of formulas for computing the probability that a point taken from 
A normal bivariate elliptical distribution with specified x, 4,0, , Sy 
shall fall within a circle cf a given radius whose center is displaced a 
given distance from the center of the distribution. The other part consists 
entirely of probability tables. The input parameters (entries) to these 
tables are given in units standardized with respect to Sy , namely: 


Ry , D/o,, , Oy | Oy . The spectra of values for these input 


Preceding Page Blank 


+ A UNI ts 


pavameters are: 
20): 0.04r £ 4.0 , in intervals of 0.2, 
R/s, :0.08R/5,£4.5 — , in intervals of 0.1 
D/s, | | 2.004 Dl, £ $00 , in intervals of 0.05. 


3. DISCUSSION 
@. THE NORMAL BIVARIATE DISTRIBUTION 


The normal (Geanwian) bivariate distribution function is given by 
the following expression: 


= ; 
f(xy )= me orl $ (32) is (<3)] ds 


where x and y are statistically independent stochastic variables (zero 
correlation), and X and ¥ are the x-mean and y-mean respectively. The -unc- 
ie f(x,y) is defined over the entire xy Blane) and dis not equal to 5 


y 
oz Sy). 


Incase Oy =O, = 0  , equation (3a-1) is transformed to: 


sage |= dolar eee] = 9059). (38-2) 


Equation (3a-2) can be converted, by means of a linear trans forma- 
tion, into one with new variables u, v, that are stcchastically independent 
and are normally distributed with mean equal to zero and with standard 
deviations oi = Sy =- cS . The transformed distribution function is 


given by: 


! a exp| - (urs v2y/0*] = h(t, vy) (3a-3) 


If x and y are not statistically independent, the distribution 
function contains the correlati-n coefficient 7 . It is given by the 
following expression: : 


F(s,y)= : wh och smc 89-9) , (4-7 ¥|\_ (30-4) 
a6 Ti-7 wim ae) ee G, ie 


The two-dimensional distribution given by equation (38-3) is obtained from 
(3a-4) by equating 7 to zero and by transforming (by a linear transformation) 
to new variables u, v. It is pertinent to observe that equation (38-3) may 
be transformed to polar coordinates to yield the distribution function 

PC p) given by the following expression: 


p(P)=£ exp (pha) for O£P<@ 


| — (3a-5) 
pl(P)=0 for ~wepco 


where is calle. the radial error: p's Wey? =(x-2)+(y-5). 


Equations (3a-5) represent the well-know Rayleigh distribution. 
The probability that a radial error 9 lies within a circle of radius R is 
ottained by integrating equation (3a-5): . 


| Rs ae Be 3 epee 
P(eer)=4, exp(-P/202)pdp =1-e —, (Ga-6) 


(4) 


It is to be noted that the distribution function (the probability density 
function), £ (x, y), in (3a-4) is constant for a constant value of the 
exponent, that is, for: 


| jend\ | Ae: Yrs y-5) . | es | 
Q(x N= ~=) -27 ae) a) =k (1-77). (3a-7) 


For a given value of the correlation coefficient 7 and & constant value of 

fe , the loci of the points (x, y) satisfying equation (3a-7) are equi- 
probability ellipses (contour ellipses). The distribution function f (x, y) 
assumes itsgreatest value at the point (5 ¥) . The equi-probability ellipses 
have the property that outside the ellipses no values of the density are 


3 2 


greater than those assumed on the interior. When 7=0, equation (3a-7) 
becomes: 


; a 2 
oa (22=z 9-7) ~ 42 3 
as (2) +($ 1 ae | (3a-8) 


Equation (3a-8) represents an ellipse for a given value of K , with its 
center at (xX, y), and with its major semi-axes parallel to the x-axis if 
Oo > Sy , and parallel to the y-axis if (25 Lg . The two principal 
semi-axes of the ellipse have magnitudes KO) and k6y respectively. The 


tables of this report will deal with ellipses represented sy equation (3a~-8). 
The probability that an error lies within this ellipse with parameter 4~ is 


2 


a 5 
pe jatoe (3a-9) 


ee a of this report will deal with ellipses represented by equation 
3a-8). 


b. OFFSET CIRCLE PROBABILITIES FOR THE NORMAL BIVARIATE ELLIPTICAL 
DISTRIBUTION. 


— Let f (x, y) be a normal bivariate distribution function with XY 
and ¥ (the means of x, y respectively) at the point M(z, gy) in the 


xy plane; and 6 ¢ 6, oo 
ea N (4,93 (593 ZFG 57 =0] . 
The probability that a point (x, y) from a normal bivariate 
elliptical distribution N[(x,y)3 (2,9); 4, ¢ #6,| with zero correlation 
falls within a region S in the «y plane is given by: 


P [44<5|= aig |e £|(=2}+ esF dx dy,( 3-1) 
o 


Let the region S be a circle @ of radius R ani the center (0, 0) at 
the offset distance D from M(Z3y)  , the center of the distribution. 


Then . 


Plewee= arb ctaf fools Carls Cet i] hea (3b-2) 


| 4i2aePer® 
Equation (3b-2) yields the probability that a point (x, y) from a normal 
bivariate elliptical distribution, N((ayu)s (259); Ox = By ;zero 


correlation | , falls within a circle £ of radius R, and ae center at (0, 0) 
is offset by a distance D from the center of the distribution: 


The following figure illustrates the configuration of the center 
M(X,¥) of the distribution relative to the origin as well as the distance 
from the center (0, 0) of the circle to point M(x, y): 


FIGURE 1 


It can be shown (see Technical Supplement) that equation (3b-2) may 
be transfurmed to: 


2 = 
Plawet]= —t- ate (| exp nil eye (ey drdy (3-3) 
Litt y eR 


This is dene by rotating the XY axes through an icangie & (= arctan Y/z) 
so that the center of the distribution is made to lie on the new X-axis (X’) 


at the point (D, 0). 


If the new xX° Y’ axes are translated through the distance D to the 
point M(x, y) and then transformed to polar coordinates, equation (3b-3) 
becomes: 


a 2 
p (or 1) cos 0+1)-2r*ppcos - 
eae 1-2 an1re.3 aca 2Piq? at dpdg (3-*) 
i P 0 ‘0 | | 


049saT 


,/2 
=—4, #2 
where r= Sy /5, and D = (z + Yy ) . (See Technical Supplement for all 
derivations) 


In order to evaluate equation (3b-4+), the integrand is expanded into 
a power series in end @ . After many transformations, substitutions, 


ca ia and regroupings, the following expressions are obtained: (See 
page 7, 


Pleeses]= Sa 


where: 


anre} 


o> 
oo 


fare a de (30-5) 
, | 
=re* mil es\, s* | CE) Ae (30-6) 
=re fe e ra [> af Gnplee "} 3b 
| Nz0 K=0 20 
" . 
eee Qn) ss\, St llan # (bjayt |) (30-7) 
= > ah j il > 
. n=O Kz0 =e 


+2 [re ‘atra][S- b] 
+ 5euj[1-€° Going Py F- St +H] 


gators als-& Gn il Sarees doa 


as 2r*pyo;=4r*t 
b=r*-1 


r= Sy/Ox | 

$= RUC = 4, (R202) (30-9) 
t= Dips?” 

A= (r?--1)c0s76 

B=2r*bdcos0 

C=2r702= aay . 


_ Note that the number of terms in each bracket of equation (3b-8) is finite, 


but the meabter of brackets ig infinite. 
a) Normal Bivariate Circular Distribution 


E If > is set equal to zero and ato in equation (3b-8) (so 
that r=1, Co, =G=6), # we obtain the offset circle probability for a 
normal bivariate circular distribution. Equation (3b-8) is now transformed 
- into: 


Play eé]= é‘{:- é — [i-e -” +#[Le “ng ] 
a a (3b-10) 


=0 

2/263 t= d/207 | 
vhere $= R*/207 ana t=TZ . This is the probability that 
a@ point (x, y) from a normal bivariate circular distribution falls within a 


circle whose center (0, 0) is offset a eae 2 from the center M (xX, jy) 
of the distribution  N[(24);(7,9);G=0j=03 zero correlaticn’) . 


(2) The Offset Distance is E to ie 


If, in equation (36-8) the carters a is set =O, but d£0, 
then D (the offset distance) = 0. Equation ey becomes: 


N-> aD 


ig | 
Plexwee]=r . secTomnce ct E a =| (3b-11) 


=0 
where $= RVAT*O? = 7/26}, 


| (3) Input Parameters 


It is desirable for many reasons that P[¢ 2 YE £] be 
dependent on j* Rb, ‘ and D/g,. . These constitute the input parameters, 
standardized with respect to Oy, for the tables; 1.e., these parameters are 
measured in units of d,. Rs (30-8), in terms of the standardized para- 
meters, becomes: 


| 4 Nooo 2 
P . £ re? * Gr) iy sane (RB) 
7 [(as)< ] Zi 1 ue eran Ki} 

=0 K=6 

(3b-12) 

an- : ies 
. eres, Ce) Dade tmious Bs 
=0 . ° 


Equation (30-12) is the infinite series representation, in terms of the 
desired standardized parameters r(=G/s,) » R/O, , andd fie » Of the right- 
hand si side of equation (3b-3). It is this equation that enables us to compute 


the probability that a point (Y,y) from a normal bivariate elliptical 
distribution with a given r(=6, raps will fall within a circle of a give 


standardized radius R/g, whose center (0,0)is at a given standardized ais 
tance tance D/e, from om the point M@, y) »the center of the distribution, where 
x > y are the x-mean and y-mean respectively. 


In order to compute Pic (Z%,Y)E€L| we enter the tables with a 
knowledge of I =0/6,) » D/5, MARA» obtain the corresponding 
probability. It is to be noted that one need not be given the coordinates 
of the center of the distribution, i.e. M(x,Y) . Only the offset distance, 
D, is sufficient for this purpose, in addition to & , oy and R (radius of 
the circle). If M& y) 18 given or determined from the distribution, D can 
be computed from the expression 


D*= x +y". (30-13) 


For ssaudie’ given: 7 (=D/6)=0.2, Rh = = 1.8, D/s, = 1.5. 

P is calculated to be equal to 0.6132. Note that were ratios are sufficient 
for evaluating the probability; one need not know the values Fd, 
D. Naturally, given the values of G,Ox ,D(or Yana G) anf, 
the entries to the tables can te made in the form of ratios. 


x ? 


Suppose ‘ we are given the problen: Compute the preoability that | 
an impact point ie & normal bivariate elliptical distributioi with 
Sy = 500 feet, Sy = 300 feet, whose center (of distribution) is displaced a 
D = 1000 feet from the aim point, should fall within a circle (about the aim 
point) having a radius of 1750 feet. In this problem r= 0.6, D/dq = 2.0, 
R /8y* 3.5. The probability is computed to be 0.9256. This means that a 
circle having a radius of 1750 feet and whose center is at the aim point 
offset 1000 feet from the center of the distribution will contain approx 
mately 92.4 of the impact points, 


Ce Eamples of the Use of the Tables Tables 


The entries (input castes to these tables aaa the ranges of 
their values were show in the summary (Section 2) to be: 


Vr (=5y/5,)2 0.0(0.2) 4.0 
Rb, — +: 0.0(0.1) 4.5 


Dio, 3: 0,00 (0.05) 5.00 . 


There is one output parameter, namely P; the tables consist of values 


of the probabilities. However, given any ““~ee parameters, one can find the 


fourth from the tables (interpolation may cimes be necessary). It is to 
be noted that in all the pre lems discusse. .clow it is assumed that the 
major and minor axes of the equiprobability ellipse are paraliel to the x 
and y axes respectively. The center of the ellipse is (x - Suppose 
we wish to consider the bombing of a target. If the aim it (or target's 


position) 1s assumed to be at the origin of the xy plane, then (x, y ) 
measures the "aim bias". 


10 


(1) # of Problems. There are several types of problens 
encountered in theory and practice, all of which stem from the fact that 
there are essentially seven parameters (not all of which are independent): 


XY 3 Ox »Fy5 R/o, , d/é, 3; P 


(2) Problem 1 Given: 6 Oy; X,y 3 R. To find P. 


We are given eriors which are governed by a normal bivariate 
elliptical distribution: Sy = 5OQ _ feet; Sy= __200_ feet; x= L,600tt. 
y = 1,209 feet. To find the probability that an error should fall 
within a circle (about the aim point located at the origin) having a radius R 
of 1450 feet. First, we compute the value of D from the formula: 
pD’?= x44 77 . Next, we compute the ratios: r=F/s5,= 200/500 BOF: 
R /o.= S00 2% 4 =; and d/6, = 2000/sag=4.09. The probability P is 


shown from the tables to be equal to 0-1299 . Tms 13 4% of the 
errors fall within a circle about the alm point and having 4 Tadius of 


1450 feet. 
(3) Problem 2 Given: 5x ys D; R. To find P. 


The impact ints of a missile are known to come from a normal 
bivariate aistrivation the mean of the distribution is displaced 3,250 


yards from the aim point; 6,= 2,009 yards; dy= 4,600 yards; to find 

the probability that an impact point will fall within a circle, about the 

aim point as center, having e radius R of 2,400 yards. In this problem, 
r=6,/6,=_1.6 3 D/e,= 32S0/lo00= 3.25 ; R/s;, = 2490/1000 = 2.4. 

‘The probability is shown from the tables to be  0./J42 . 


(4) Problem 3 Given: 5 Sy; dD; Pp. To find R. 


. What is the radius, R, of a circle, about the aim point as 
center, that contains % of the ct points of an artillery shell? 
The impact points are governed by a normal bivariate distribution: 


Ox= 800 feet; Sy= 3,040 feet; the mean of the distribution is 
displaced £560 feet from the aim point. First, we compute the ratios: 


r= dy /S, end D/6x . Next, from the tables, w determine the 


value of R /& 4-0 corresponding to the values of P 9 D/ox » and 
fT. Finally, we obtain the value of R = 3,200 feet from the relation: 


R= (R/,) & 


(Note: f= 3040/g00=3.8; 2/6, =1560/800=1.95 ; Pz0.60). 
aBe 


(5) Problem 4 Given: 6,6,3 R; P. To find D. 


The (two-dimensional) errors of a target's position are governed 
by a normal bivariate distribution, with: oy * S00 _ feet; Oy = 400 _ 
feet. 50 % of the errors lie within a circle about the true position as 
center, having a radius of 550 feet. How far is the mean of the distri- 
bution displaced from the center of the circle? Pirst, we compute the. 
ratios: r= y /Sx ; R/dx + Next, from the tables, we determine the 
value of D/6, = 0.40 _ corresponding to the values of P, R/S, » and r. 
Finally, we obtain the value of the displacement distance, D, from the 


relation: D=(D 6, . (Note: £ = 400/500= 0.8; R/6,= 550/s0pz1.1 
= §P/sx) . ( e i A 1, “ Bh, 
(6) Problem 5 Given: Pe R; D; Sy . To find Sy. 


We are given target position errors which are governzd by a 
normal bivariate elliptical distribution. We are also given the probability, 
Pe bs that an error should fall within a circle (about the aim point 
as center ving a radius, R, of 2,850 feet and whose center is displaced 
@ distance D =- LL25 feet from the mean of the distribution. If, in 
addition, GO, is given equal to 750 feet, whet is the value of dy ” 


First, we compute the ratios: R /Sy. 5 D/ & . Next, from the tables, we 
determine the value of ff (= by /5x%) * 2.§ _ corresponding to the value: 
of Py, R/é,y and D/s,. Finally, we obtain the value of Oy from the 


relation: 
By =T6, = 2,100 fee, 
(Note: R/6, = 2850/750= 3.8; D/y = 1125/750 = 1.50) 
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Se TECHNICAL SUPPLEMENT | 


. In this section we shall show how to obtain equations (3b-5), ( scsi 


ae és lTot FA+-07| p dp do (5-1) 
a eel fo Ber on 
=re* ) {aa ey s'](e "ye nee tela} (5-3) 


K=0 


where the parameters A, B, C, a, b, r, 8, t are defined in (36-9). We see 
from equation (3b-1) that the general expression for Ply) 5] 
is given by 


| oka, pe! 
Plansl= te | Jeol (8 a Ce I dx dy os 
. 


where the correlation coefficient is zero, the region is S, and X and y 
are the means of x and y respectively. Suppose we wish to determine the 
probab! .ity that a point taken from a normal bivariate elliptical distribu- 


tion N[ Cayrs Uy) § (6, Sy) shall fall within a circle of given radius R 
whose center is displaced @& given distance D from the certer of the distri- 
bution, ™ (%¥)- Equation (5-4) becomes 


Pibondless| oo alt h+Cylfaa. on 


e+ y"£ Rt 


13 


In order to evaluate the double integral in (5-5) we shall rotate the 


coordinate axes through an angle &f= tan '( g/x ) and translate to 
Mm (X,¥) through the distance D (See Figure I). We obtain | 


: ae ak 1|/x-D\ | 
| Plaweel= =A J exp {252 + @) i dxdy. (5-6) 
x44y%R? 
Next, we shell transform (x, y) to polar coordinates: 
x=pcos6 | t 


y= p sind (5-7) 


eee 


0<p<R (5-8) 


Then 


ddy=|J(2)| dpde=pdpde. 2) 


Now set | . rr by 3 (5-10) 


Equation (5-6) becomes, after simplification, 


Diao? (27 -R (5-11) 
Ple,»ee]=<— | \ e pdpde | 


ik 


vhere 
| y= LA+1)- PB 


A= (r?-1)cos7@ 
B= 2r*p cosé (5-12) 
C= 2r%s2 = 26," 


To evaluate the double integral ir (5-11) we use the infinite series 


1 | | | 
5 i uf | | (5-13) 
ae me + we | 


We employ the following definite integrals: 


TY an ar an+ 
- 3 ae << = 
NN f cose de = { cos Odo =-<- ={ Cos bdo=o 


2sr a 
i cos*6 de=mr , J cosnede =0 7 n>od 
(+) 


an — (5-14) 
| an 305-6 Teoe Ca-)) > 7 
Soret rereeenemeee 


—~ 3:5-7-.. (n-1) 
? 2"-in| Tr 


Substituting ed) in (5-11), we get 


af cy par -R 
ue 
Playet]= 2 


Soraz | \(- fee Be fel dapde «can 
o 
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To evaluate the double integral, we integrate term by term, resorting to the 
binomial expansion of U : 


n n+t TN 
f= “(Pl +1)- B] a (5-16) 
12 ae ie tants 
pues e (7)0 ase BS 
(=0 


=a “> ye ) OCP) A tele mV an 


By carrying out the integrations and applying ieany we obtain the infinite 
series (see page 17) 
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Ployec] 


Re 2. 2b . 3b* _3.3ba?_ 3a*, 3a 
6 G [getaet ai} waicd 31C8 Hues 


ab P 35d) _ 6a" _ 3.1zba® 
3!¢3 seat ree H6-310> Her gH 


3.5-6bu*  3.5.5ba* 35a 3-50 


~ FEM Ce” HE5ICS | H5ICS  46-6!C8 ra 


RYT 2, ub , 3-667, 3.5.4b? 3.5.76" 
uice  wict wwict 4.6-4ICt #6.8.HICH 


100? 3-3.10ba* 3-3-5.10bK* _3.5.7.10b4" 


5icS 515 46-5108 — 4.6-B-51C5 


3.1504 35-15-2ba*  3.5.7-15b*) 
a #-6ic8 4.6-6!1C° a 4.6-8-61¢c& 


(5-19) 


In equation (5-18), let us put 


a= 3 ¢ (5-20) 
Then, from (5-19) and (5-20) 


* 7 2 
x 
Tus aort-ac. ' (5-21) 


Referring to equation (5-18), we see that it involves even powers of 

ow: X27, at ,..., so that | | | 
aYC=a; a4c*=al; x4c*=a'/c*;..-;07"C"= aE >m) 

Substituting (5-21) in (5-18) one obtains an expression for P[(x%yéL 

involving D/s, , r, Sy, »,andR. Tms, lay ] 


Pl yee] 
e meet eee -£,] 
arai) 2 # [cC "e ~21¢ 
+E 2.,2b 4 3b _ 3.3ba_ 30, 3a? 
6 [2c * ic® * hatc® — Gace 30,4 208 
R 


sv 2 3 7 
2, 3b | 3-36, 3-56 6a _ 3.12ba | 
E> 3+ y3Is THEO Hid HuICS (5-22) 
_ 3:5-66@, 3.5.5ba*, 3.50" 3.5 
4.6.4) C3 4.6.51¢C3 4.5! Cc? 4.6, 6! Cc 


+ Ho , 3.66, 3-5.4b>, 3.5.76" 
10 [Mice wet uict 4.6.44 4,6.8.41¢C% 


_ 10a 3:3-loba _ 3.3-5-10ba 3-5. T-loba. 
5ic* 4.5ic# 4.6-51C% 4.6.8-51¢4% 


3-150, 3-5-15-2ba* | 3.5.7-15b0"|_ 
4.61C% 4.6.6IC+  4.6.8.61C% 


In (5-22) let us make the substitution 
me y/ 
aw a a (5-23) 
80 that R'/C=Cs* Ric*= Ce Gk JC" '=CS* We then obtain the equation 
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Plyec] 
_ ce* 


are 


$*[, al, #i2 3b* _3.3be _3a 
# [+ $-|+ 6 z+ a1 ar 2! R+#- eat 3! rr 


-#[2 +264 2H" ssb*, 35b? 6a 3.izba 3.566% 
4.3! + e603! Gt bl 4.6! 


1 3:5-5ba®, 350° _ 3.502 
+651 * ST #66! 


(5-2h) 
5 > % 
s Mb , 3.667 35-40", 3. 3.5.76 10a 3-3-10ba 
tit ae +a rere yal * heat Wu! 4684! ST yg 


_ 3-3-5-10b0 _ 3.5-71oba , 3 15. a 3-5-15-2ba* 
4.6.5! 4.6-86-5! 4.6! 4-6-6! 


4g 2:57:18 BO" 3-5-10?_3.5-7-7b08]_ 
ee 468-6! 4.6.1! 468-7! : 


In (5-24) we shall combine terms in such a way as to lead to a pattern or 
law of formation. We shall group terms within the brackets and add vertical- 
ly. Adding the first colum of terms we obtain 


2 3 rT § 
Ge ers Pol 8s ts 8 : 
sels Ua ads ee | = 
It is well kmow that 
e= “itar-artgr-gart: (5-26) 


Thus, we see that (5-25) is precisel: 


-e°= [ -é*][1] : | (5-27) 


Again, in (5-24), if we add vertically the terms involving the first powers 


of a and b, we obtain » 


-3¢ wm 066 (5-28) 


$-9-£0-8). +3 


= --£( -F)+e s’ 


Collecting and simplifying the expressions in (5-28), we obtain 


5 & 5 
1 s* 28°, 38" us : 
In (5-29), the expression in the second parentheses is precisely 
=—§ 
t-e (s+1)| . 
Hence the terms involving the first powers of a and b> add up to 
(5-30) 


oz -4[t-eF(1+5)][2—-b] - 


Again, in (5-24), if we add vertically the quadratic terms in a and b, t.e., 
the terms in a@, ba, and b "y we get 


aeeenennetntt tee stmediiantien oinioel 


Sum (quadratic in a and b) 


— $3] 3a*_ 33ba_, 3b* 
~ 6 Laut mt” Gal 


_ SHlasa* 3.12ba , 3.3% 

81451 ep 

F | 
_ S°| 3150" 3-3-40ba | | 
+§ [sae rs | : | (5-32) 
_ | 8-352" _ 3-4.15ba | 3-10b* 

12 | 4-7 +6) cri 


7 
$'| 3-700" 3.21. 3:21-5ba , 3-15b 
% [ze “en + se 


By collecting and simplifying the expressions in (5-31) ve obtain 


2113] s* sha sF_ 4 56457 [a2 sha, bt] (5-32) 
sun ofA ora 10 S!iz aie wa) tar]: 


: (5-33) 


Hence 


Sum (quadratic ine a b) 


= 213 i-e° 1+ + a - she 
2-4 i i ah RE es 
2-4 E e alle = oT +#] o : ss (5-35) 


Continuing in this manner, we obtein 


Sum (cubic in a and eae ety” *l[¢- Shed +S fl (5=36) 
2 Cy 
Sum (quartic in a and vetiee re royale ~ he sie Smelt 37) 
5 K § % 3 
541357442 a? abe", 28be? astia® 
i agar “Hasle Zi ‘vals at gl ~ iste 


Let us derive the expression for the general or nt term. wWe shall use 
equation (5-38). The derivation, however, is perfectly general. The first 
term of the right-hand side of (5-38) is 


(535-29 0” : 
a "246810 9 a 
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ieeestenenhiatines Sainte 


The denominator in (5-39) can be shown to be equal to Z-5! , so that 


(5) 9°51.1.3.5.7.4 _ 101 (2-5)! | | 
Ty. ee eas ee ‘ (5-40) 
210, 5) 210.5; 28s! | 


Proceeding in this manner, it can be show that 


(m)_ _Q@n)l — 1-3-5... (2n-1) _ : “4 
agp =Sar, = mer | eae (5-41) 


Next, in (5-38), consider the 2°¢ tern, — , in brackets: 
. § 
5) -$ K 
Be[-e re] (542) 
The general term , Ts is obvious: 


n) a | 
=O 


Finally, consider the gra expression in brackets on the right-hand side of 
(5-38): 


()_ a qbab 2eba asba?, isb’a ob : 
Zs Toi” “qa 7 6! - =n ter Br ; oe 


is aenreatare as (S-hh) are: 
en! 5 (RSD! (25-2)!5 Q-5-3!,(R5—¥)!; (25-5)! 
The numerators in (5-44) are 


crja’, »(-1) Cc 5-1 " oot. i (* ae ft, 1 He syne 5-3, 


C1)'(2 5 )plars (-1yb° . 


Tms, each term in (5-4) may be written: 


i (* we 
bén 5 -$ 
(-1) oa : . 42041, 195 (5-45) 


Equation (5-38) may be written in compact form as follows: 


| + (3 ih ‘4 : 
Sum(Quintic in a and vegan fall ye |. (5-46) 


@ Dl 


We can now write down the general expression for the n-tic in a and b. It is 


ie eS ot CE :(™ CE) ed (5-147) 


(2n-4)! 


ae 


intent een a ed 


Sia; pigté = *(b/a)t 
en | CO) (an-$)! _ G@n- 


a Qn-j-f! A feen-pl ° 


we may write (5-47) in the alternate form 


2a, Ars | 
yn C F- Hain “pl (5-48) 


Combining equations (5-41), (5-43), (5-47), and (5-48), we obtain the required 
probability Pla ye 7AR 


Plomedlesdcl four {-+[s2h+(2i]} aa 


Ex “y"¢R? 


-Bhe (2 ee 
=te | lor {- a [(r%1) cos a dp de 


=e) a ofales “tea 


“\bere, in (5-49) and(5-50), 


saz (8) ;r= 2 
¢ +eY 5 tent 1-3-5.. ae 


2™n! (5-51) 


a= ar*(D/s,). 
b=rt- 1 


Note that the numter of terms in each bracket of (5-49) and (5-50) is finite. 


EL a ses a ee 
Playee] 

=ré ‘{i-¢ in 

i [1-€*G+9][F - 5] 


+H eet +S e- Bey ¥ | | (5-52) 
+;£4(1- eo slg- Sba* fe #] | 
al Rae ee ah 

lve gas Bae Sea g 

+: + Blt € ‘yt S-a i + CMa” 7 ae), 


26 


Let us make the following substitutions: 


A= R/, | 
Lu = D/s, : (5-53) 


Then s, t, and a in (5-51) become 
$= X/2r* 
t= | _ (5-54) 


ax arth 


‘Substituting (5-54) in (5-49), we obtain: 


PLove Al 


aS (ge es i Salat a eh. roy 


N=0 K=0 


Equation (5-55) is identical to equation (3b-12), and involves the standard- 


ized parameters: lr (=4,/s,) - R/6, , ‘ —* . Let us now consider the 
following two special cases: 


Ge az 0, b= O~. 
Setting b = 0 in (5-51), we see that P= 1 : 6x =6y=6. 
Also, $=R'/26%, t=D/267, a=2(p/s)*= 2°/267)= 2t. 
y, 2m)! | am NN +t" Thus, for b=0, 
zn! Qni 24% n! n! | 
we have: 
6x = 6y=6 
$s = WRG 
ee ees 
t= vie eax) 
aS :2 
Qn)! a® — t” 
2™n! (an) a i 
Substituting (Sa-1) in (5-49) @r (5-52), we get. 
N»2ao n , 
-t n 7 Ki 
Pl@weel=e if rey sf ? (5a-2) 
M=0 are : 


From equation (5-6), we see that when r =1 (i.e. when 6,=dy= co), 


Ge 2), ¢) )=4 [o-py+ y7] : (5a-3) 


The right-hand side of (Sa-3) equal to a constant represents a circle. | Thus, 
when b » 0, the distribution is circular. Hence. P (x,y) é £] 
_represents the probability that a point (x, y) from normal bivariate 
 Gareular distribution, N[(xy)3 9) 360) falls within a circle 
£ of radius R whose center (0, 0) is offset a distance D from the center of 
the distribution. ae ee oe . to . 


b. a= 0, b#0. From equation (5-51) we see that vhen a = 0, D = 0, 
since r0. Thus, the offset distance is zero, and the mean (X,J) of the 
normal bivariate elliptical distribution coincides with the center (0, 0) of 
the circular region £& . . . 


Substituting a = 0 in (5-51) ve obtain 


D=0 | 
t=1(D/5,) =0 } 3 (50-2) 


Putting the relations (Sp1) into (5-49) or (5-52), ve obtain 
| a = (-1)" (an)! (r21) . -$ 5K 
Plowee]=r | oe e =| : (5b-2) 
m0 oe 


Thus, when a = 0, we see that Play £| represents the 
' probability that a point (x, y) from a normal bivariate elliptical distribu- 


tion, N[G.y5 (0,0); (5, ,6y) | falls within a circle £ of radius R whose 
center (0, 0) coincides with the center of the distribution. 
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